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Abstract 

Free induction decay (FID) measured by nuclear magnetic resonance (NMR) in a polycrystalline 
solid is the isotropic average of the FIDs for individual single crystallites. It has been recently pro- 
posed theoretically and verified experimentally that the long-time behavior of single-crystal FIDs 
has the universal form of exponentially decaying sinusoidal oscillations. Polycrystalline averaging 
complicates the situation theoretically, while the available experimental evidence is also ambigu- 
ous. Exponentially decaying sinusoidal oscillations have been observed for 129 Xe in polycrystalline 
solid xenon but not for 19 F in the powder of CaF2- In this paper, we present the first principles 
FID calculations for the powders of both CaF2 and solid xenon. In both cases, the asymptotic 
long-time behavior has the expected form of exponentially decaying sinusoidal oscillations, which is 
determined by the single crystallite FID with the slowest exponential decay. However, this behavior 
appears only at rather small values of the signal that have not yet been measured in experiments. 
At intermediate times accessible experimentally, a polycrystalline FID depends on the distribution 
of the exponential decay constants and oscillation frequencies for single crystallite FIDs. In CaF2, 
these parameters are relatively broadly distributed, and as a result, the sinusoidal long-time os- 
cillations become somewhat washed out. In contrast, the single crystallite parameters are more 
clustered in solid xenon, and, as a result, the experimentally observable range is characterized 
by well-defined oscillation frequency and exponential decay constant even though both of these 
parameters do not represent the true long-time behavior. The above difference of the intermediate 
FID behavior originates from the difference of the crystal structures of solid xenon and CaF2- 
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INTRODUCTION 



First principles calculations of the free induction decay (FID) measured by nuclear mag- 
netic resonance (NMR) in solids is a long-standing theoretical problem [1] still lacking a 
controllable solution 1]. The most challenging aspect of this problem is the prediction of 
the long-time behavior of the FIDs. Recently some progress in this direction was made on 



the basis of the notion of microscopic chaos [2 



-|4|. Namely, it was predicted that the generic 



long-time behavior of FIDs in single crystals has the character of exponential decay with 
or without sinusoidal oscillations. In the most common case of magnetic dipolar interac- 
tion between nuclear spins, the oscillatory regime is realized, and hence, the long-time FID 
behavior can be parameterized as 



F(t) = Ae~^cos(ujt + 



(1) 



where A, 7, u and <j> are some constants whose values were not predicted. It was only 
estimated |3| that, generically, the values of 7 and u fall on the timescale of the spin-spin 
interactions often referred to as Ti- It was also estimated that the long-time behavior (pQ) 
becomes dominant after a time on the order of several times T2 from the beginning of the 
FID. The above predictions agree with the experimental 5|47| and numerical[8|, |9j results for 
quantum and classical spin systems. 

The situation becomes somewhat more involved theoretically for polycrystalline samples 
or crystal powders. Different orientations of single crystallites in polycrystals/powders with 
respect to an external magnetic field imply different microscopic Hamiltonians, and hence 
different values of 7 and u, which in turn leads to the additional averaging over the oscillation 
frequencies. At sufficiently long times, the crystallites exhibiting the smallest value of 7 
should start dominating the overall response, and, therefore, the well-defined frequency 
of these crystallites should also control the overall decay. We call the latter regime the 
asymptotic long-time behavior. It is to be distinguished from the intermediate behavior, 
which we define as the regime, when the individual crystallites have reached their respective 
long-time regimes but the asymptotic polycrystalline long-time behavior is not yet reached. 
The challenge here is to understand how long the above transition to the asymptotic behavior 
takes, and what the intermediate behavior looks like. It is, in particular, possible that the 
intermediate behavior exhibits a tentative "washing out" of the FID beats. 
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On the experimental side, the available facts about the long-time FID behavior in poly- 
crystals/powders do not reveal a consistent picture. On the one hand, no well-defined 
long-time beats of form ([T]) have been observed in the CaF 2 powder (within the range lim- 
ited by the experimental signal-to- noise ratio) jo[ [lo| . On the other hand, in hyperpolarized 
solid xenon, which is supposedly p olycryst alline , the experiments reveal well-defined beats 
of form (pQ) appearing rather quickly |6l. Ill|. 

In the latter case, the situation is complicated by the fact that hyperpolarized solid xenon 



is prepared in convection cells 



12| by first optically polarizing xenon gas 13) and then rapidly 



cooling it into a liquid phase and subsequently quenching the liquid into the solid phase. As 
a result, an uncertainty remains about the proper thermalization of the resulting solid. In 
addition, the formation of crystal structure in solid xenon is controlled by the relatively weak 
van der Waals interaction, which is known to allow significant residual atomic motion 14 1 



that further complicates the theoretical analysis. A related unclear issue is the strength of 
the exchange coupling between xenon nuclei. 

In this paper, we assume that hyperpolarized solid xenon samples investigated in Refs 
11] can be described as polycrystalline fee lattices of immobile nuclear spins coupled by 



magnetic dipole interaction. We perform the first principles calculations of Xe FID on 
the basis of the approximation procedure introduced in Refs.j^, 15]. We also perform the 



first principles 19 F FID calculation for the powder of CaF 2 , where 19 F nuclei form a simple 
cubic lattice. Our goal is to verify whether the above calculations are sufficient to exp lain 



why the well-defined beats of form (JTJ) were seen in polycrystalline solid xenon [6|, lU[ but 
not in CaF 2 powder f], [lfj ■ 



THEORETICAL APPROXIMATION SCHEME 

We will use the approximation scheme for FID calculations that was introduced in Ref . fl5 1 
with small modifications added in Ref. Qj. This scheme is quite similar to the one introduced 
earlier in Ref. 17]. Alternative attempts to calculate powder FIDs were made in Refs. If], 18]. 



The approximation technique of Ref. [15j results in a very accurate description of the ex- 
tended initial behavior of single crystal FIDs in CaF 2 . It also leads to the long-time behavior 
of form ([!]), but with constants noticeably different from those observed experimentally (see 
below). As explained in Ref.[3|, an accurate prediction of the parameters in Eq.flT]) is not ex- 
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pected here due to the oversimplified nature of the approximation. We are, however, mainly 
interested in the qualitative question of the difference the solid xenon and the CaF2 powders 
posed at the end of the preceding section. Answering this question presumably depends on 
the qualitative differences in the distributions of 7 and u for different orientations of single 
crystallites in the external magnetic field. The approximations used should, therefore, be 
adequate for detecting such differences, if they exist. 

In CaF2, 19 F nuclei are characterized by spin 1/2, gyromagnetic ratio 7 g = 25166.2 s _1 G _1 
and abundance v = 1.0. These nuclei form simple cubic lattice with period d = 2.723 A 
(at 293 K). For solid xenon, we perform the calculation for the fee lattice with the nearest 



neighbor distance d = 4.4 A and abundance v = 0.86 o 



129 Xe nuclei. This abundance is 



111 ] . Other nuclear isotopes present 



representative of the sample most studied in Refs.jfj, 
in this xenon sample are assumed to be non-magnetic. (Here, in particular, we neglect the 
contribution of the magnetic isotope 131 Xe, which has spin 3/2 with a smaller gyromagnetic 
ratio. Its abundance is 2 per cent in the sample analyzed.) The 129 Xe nuclei have spin 1/2 
with gyromagnetic ratio 7 g = 7452.11 s~ x G _ . 

We obtain the powder FID as the average over large number of single crystallite FIDs. 
The orientation of each crystallite in the external magnetic field is selected randomly. 

For each crystallite, we calculate the FID as the infinite temperature correlation 
function [l] 

F W = — " T TX 2 X ~ ( 2 ) 

for the microscopic Hamiltonian of the truncated magnetic dipole interaction in the Larmor 
rotating reference frame: 



mn 

m<n 



p p _ -(I X I X + p p) 



(3) 



where ra and n are the lattice site indices, l & m is the operator of the 5th (x, y, or z) component 
of the rath nuclear spin 1/2 with the z-axis chosen along the direction of the external static 
magnetic field, and J mn are the coupling constants given by 

7g 2 ft 2 (l-3cos 2 fl mw ) 

— l_ _ r 13 " ' ' 

1 1 m 1 n | 

Here, r m is the position vector of the rath nucleus, and 9 mn is the angle between vector 
(r m — r n ) and the z-axis. 



Extending the approximation scheme of Refs. 0, Q to the case of isotopic abundance 
v < 1, we obtain the FID function F(t) as the numerical solution of the following integral 
equation: 



F(t) = g(t) + a \ F(t - t')^^dt\ 



where 



g(t) = exp 
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(see also |l9j). 
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The initial condition for Eq.((5]) is F(0) = 1. Parameter a given by Eq.((7j) does not depend 
on r\. Therefore, one can first set r\ = 1, then calculate a and finally use Eq.flH]) to calculate 
the actual value of rj. 



RESULTS AND DISCUSSION 



In order to illustrate the performance of the above approximation scheme, we show in 
Fig. [1] the results of the calculations of the initial and the long-time behavior of single 
crystal FIDs in CaF 2 for three directions of the external magnetic field. While the linear 
plots (insets) in each panel of Fig. [TJ illustrate that the overall agreement of the theoretical 
and the experimental curves is very good, the semilog plots (main panels) amplify the 
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[100] 


[110] 


[111] 


7, fjs 1 


UJ, [IS 1 


7, /xs 1 


UJ, /US 1 


7, fjs 1 


UJ, [IS 1 


exp. 


0.054 


0.156 


0.042 


0.103 


0.029 


0.066 


th. 


0.071 


0.143 


0.054 


0.085 


0.038 


0.056 



TABLE I: Table summarizing the experimental and theoretical values of parameters 7 and u> for 
the single crystal CaF2. The parameters are obtained by fitting the long-time tails of the FIDs 
presented in Fig. [Qby Eq.([T]). 





CaF 2 


solid Xe 


7, ^s 1 


UJ, [IS 1 


7, ms 1 


uj, ms 1 


theory (asymptotic) 


0.040 


0.055 


1.35 


2.66 


theory (intermediate) 


0.050 


0.053 


1.55 


2.48 


experiment 






1.25 


2.06 



TABLE II: Table summarizing the long-time parameters 7 and uj corresponding to the theoretical 
calculations and experiments for CaF 2 powder and polycrystalline solid xenon. The numbers in 
the "theory (asymptotic)" row are obtained from fitting Eq.([T]) to the true theoretical long-time 
behavior as exhibited in Figs.[2jb) and[3jb). The numbers in the "theory (intermediate)" row are 
obtained from fitting the intermediate behavior of the theoretical FIDs in Figs. [2jc) and Etc) by 
Eq.(H|). The experimental numbers for solid xenon are cited from Ref. Q]. 

discrepancy in the long-time tails. The comparison of the theoretical and the experimental 
values of 7 and u for the long-time fits of form (OQ) is presented in Table [H It indicates a 
typical discrepancy of about 20 percent. 

The calculated FIDs for the CaF 2 powder and for polycrystalline solid xenon are shown 
in Figs. |2] and [3j These FIDs were obtained as the average over 1000 randomly oriented 
single crystallites. 

Let us first examine the discrepancies between the theoretical calculations and the experi- 
mental curves. For the CaF 2 powder, the discrepancies appear starting from the intermediate 
section of the FID. We believe that these discrepancies are due to the limitations of the the- 
oretical approximation scheme based on Eq.©. On the other hand, the discrepancy for the 
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200 



FIG. 1: (Color online) Single crystal FIDs for CaF2- The directions of the external magnetic fields 
are indicated above the plots. Solid red lines represent the experimental results of Engelsberg and 
Lowe[5]. Dashed blue lines represent the result of theoretical calculations based on Eq.([5]). Main 

panels contain semi-logarithmic plots. Insets — linear plots. 
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FIG. 2: (Color online) F FID for CaF2 powder. The solid red line represents the experiment of 
Ref. 6]. The solid magenta line represents the earlier experiment of Barnaal and Lowe[lo|. (The 
data points are actually extracted from Ref.fli]].) The dashed blue line represents the theoretical 
calculation described in the text, (a) Linear plot, (b) Semi-logarithmic plot with the thin green 
line representing the fit of formal]) for the asymptotic long-time behavior of the theoretical FID. 
(c) Semi- logarithmic plot with an attempt to fit the intermediate FID behavior by formula ([I]). 
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FIG. 3: (Color online) 129 Xe FID for polycrystalline xenon. The solid red line represents the 
experiment of Ref.Jg], the dashed blue line is the theoretical calculation described in the text, (a) 
Linear plot, (b) Semi-logarithmic plot with the thin green line representing the fit of formflT]) to the 
asymptotic long-time behavior of the theoretical FID. (c) Semi-logarithmic plot with an attempt 
to fit the intermediate FID behavior by formula (JTJ). 
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solid xenon appears from the very beginning of the FID. It is related to the fact that the 
theoretical and the experimental values of the second moment M 2 are different from each 
other (2.64 ms -2 and 1.6 ms -2 , respectively). Since the theoretical value of M 2 is the in- 
put rather than the output parameter for the theoretical approximation scheme, the above 
discrepancy indicates the inadequacy of our initial assumptions about either the form or 
the parameters of the Hamiltonian ([3]). It may be related to the insufficient thermalisation 
and/or atomic motions in the quenched solid xenon samples j^]. Leaving this discrepancy to 
be investigated in a later experimental work, below we focus on the outcome of the theoret- 
ical calculation and examine the differences between the long-time FID behavior for CaF 2 
powder and polycrystalline solid xenon. 

Figures [2^b) and^b) include fits of the true theoretical long-time behavior to the asymp- 
totic formula flTJ), while Figs.^c) and[3^c) attempt to fit the middle section of the theoretical 
FIDs with Eq.flTJ). In CaF 2 , the intermediate FID behavior is not well described by Eq.flTJ). 
At the same time, the asymptotic long-time behavior becomes pronounced relatively quickly 
- after about three beats. On the contrary, the intermediate behavior of the solid xenon 
FID is well described by Eq.flT]), which covers about 6 beats and 5 orders of magnitude, while 
the asymptotic long-time behavior emerges only at relatively late times and small values of 
FID. 

It is expected that the behavior of the middle section of the FIDs is controlled by the 
typical single crystallite values of 7, while the true long-time behavior is controlled by the 
crystallites with the smallest value of 7. In order to clarify this issue further, we present in 
Figs. E|a) and E(a) the theoretical values of the long-time parameters 7 and u for the single 
crystallites included in the powder average, while Figs. Efb) and[5^b) show the histogram of 
the resulting points. 

One can now appreciate the qualitative difference between the powder of simple cubic 
crystallites and the powder of fee crystallites. The long-time parameters 7 and u for the 
simple cubic lattice are much broader and increase or decrease roughly proportionally to 
each other. Therefore, the typical values of 7 and u are sufficiently different from those 
representing the asymptotic decay. The poor performance of the middle section fit is in large 
part due to the larger difference of frequency u) between a typical value and the asymptotic 
long-time value. 

On the contrary, all possible values of 7 and w are more clustered for the fee polycrystal 
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FIG. 4: (Color online) (a) Theoretical values of parameters 7 and u for 1000 randomly chosen single 
crystal orientations of CaF2- The blue shape appearing in the plot consists of 1000 points. Each 
point represents a pair of values (7,w) for one single crystal orientation, (b) Histogram iV/j(7,w) 
of all sampled points in (a). 
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FIG. 5: (Color online) (a) Theoretical values of parameters 7 and u for 1000 randomly chosen single 
crystal orientations of solid xenon. The blue shape appearing in the plot consists of 1000 points. 
Each point represents a pair of values (7,^) for one single crystal orientation, (b) Histogram 
^(7,0;) of all sampled points in (a). 
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and do not exhibit much of a systematic dependence on each other. As a result, the typical 
value of 7 is rather close to the true long-time value. This explains why the fit (JT]) to the 
intermediate FID behavior works so well over an extended time interval. 

The more clustered behavior of parameters 7 and u for the fee powder was, in fact, 
expected. The differences in 7 and uo originate from the differences in the truncated Hamil- 
tonians for different orientations of the magnetic field with respect to single crystallites. The 
orientation-dependent differences are expected to be smaller for the fee lattice, because the 
fee lattice is in a sense more isotropic: each lattice site has 12 nearest neighbors as opposed 
to 6 nearest neighbors in the case of simple cubic lattice. The 12-neighbor environment is 
obviously more isotropic than 6-neighbor environment. The higher sensitivity of the simple 
cubic lattice to different orientations of the magnetic field can be illustrated by the example 
of the magnetic field oriented along [111] crystal direction, in which case the coupling con- 
stants 01]) to all six nearest neighbors become equal to zero — the so called "magic angle" 
condition. 

In principle, the polycrystal/powder average also depends on the distribution of param- 
eters A and <j) i n Eq.flT]), but we found that the parameter A has comparable values for all 
orientations and that its distribution does not add any new qualitative insight to the above 
discussion. Likewise, we were not able to find any particularly important aspect associated 
with the distribution of (J), apart from the observation that it makes the frequency u of the 
intermediate section fit for CaF2 powder smaller than the minimal value of u for individual 
single crystallites. 

CONCLUSIONS 

We have presented first principles FID calculations for the powder of CaF2 and for poly- 
crystalline solid xenon. The long-time FID decay for powders/polycrystals is the superpo- 
sition of the long time decays for individual single crystallites. The typical single crystallite 
values of the long-time parameter 7 control the middle section of the resulting FIDs, whereas 
the true long-time behavior is controlled by single crystallites with the smallest value of 7. 
We have found that the single crystallite parameters 7 and to are rather broadly distributed 
CaF2, and as a result, the intermediate section beats become washed out and relatively 
quickly evolve to the asymptotic long-time behavior. Such a behavior might be observable 
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in the future CaF2 powder experiments with improved signal-to-noise ratio. On the con- 
trary, in the case of solid xenon, the single crystallite values are more clustered, and as a 
result the middle section is characterized by well defined beat frequency and exponential 
decay constant over several orders of magnitude, while the true long-time behavior appears 
only at relatively later times. We explain the above clustering of parameters 7 and ui by the 
more isotropic character of the fee lattice in comparison with the simple cubic lattice. Our 
findings suggests that the experiments conducted so far in solid xenon have been able to 
access only the intermediate section of the powder/polycrystalline FIDs, and hence observed 
the well-defined behavior ([I]). 

It is clear that, although observing well-defined behavior ([T]) in the intermediate FID sec- 
tion requires suitable crystal structures, such a behavior would be extremely unlikely, if the 
long-time behavior of single crystallites were different from ([1]). Therefore, the experiments 
accessing the intermediate section of FIDs in polycrystalline fee solids are appropriate to 
test the theoretical long-time predictions 3|, |4j originally made mostly for single crystals. As 
discussed in Ref.J^], the same conclusion is likely true for solids with disordered arrange- 
ments of magnetic nuclear sites, but further experimental and theoretical investigation of 
this situation is necessary. 
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